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Abstract 

Viscous fingering occurs in the flow of two immiscible, viscous flu- 
ids between the plates of a Hele-Shaw cell. Due to pressure gradients or 
gravity, the initially planar interface separating the two fluids under- 
goes a Saffman- Taylor instability and develops finger-like structures. 
When one of the fluids is a ferrofluid and a perpendicular magnetic 
field is applied, the labyrinthine instability supplements the usual vis- 
cous fingering instability, resulting in visually striking, complex pat- 
terns. We consider this problem in a rectangular flow geometry using a 
perturbative mode-coupling analysis. We deduce two general results: 
viscosity contrast between the fluids drives interface asymmetry, with 
no contribution from magnetic forces; magnetic repulsion within the 
ferrofluid generates finger tip-splitting, which is absent in the rectan- 
gular geometry for ordinary fluids. 



PACS number(s): 47.20. Gv, 47.20.Ma, 47.20.ky, 03.40. Gc 



1 Introduction 



The Saff man- Taylor problem |T] , in which two immiscible, viscous fluids move 
in a narrow space between the parallel plates of a Hele-Shaw cell, is a widely 
studied example of hydrodynamic pattern formation where interfacial insta- 
bilities grow and evolve J2J. The initially flat interface separating the two 
fluids is destabilized by either a pressure gradient advancing the less vis- 
cous fluid against the more viscous one, or by gravity coupling to a density 
difference between the fluids. 

Ferrofluids, colloidal suspensions of microscopic permanent magnets, re- 
spond paramagnetically to applied fields ||. Because they are liquids, they 
flow in response to magnetic forces. Ferrofluids confined within Hele-Shaw 
cells exhibit interesting interfacial instabilities. One of the most beautiful, 
the labyrinthine instability, occurs when a magnetic field is applied perpen- 
dicular to the Hele-Shaw cell. Elements of magnetized liquid repel each other, 
creating highly branched, intricately fingered structures. 

Recent experiments ||] examine the Saffman- Taylor instability, when one 
of the two fluids is a ferrofluid, in the presence of a perpendicular magnetic 
field. The resulting interfacial patterns, in a rectangular Hele-Shaw cell, are 
an intriguing superposition of familiar forms from ordinary viscous fingering 
and labyrinthine patterns. Two immediately striking features of the patterns: 
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The pattern of low viscosity fluid penetrating into high viscosity fluid is 
totally unlike the pattern of high viscosity fluid penetrating into low viscosity 
fluid; Finger splitting is prevalent, while in general it is completely absent in 
zero external magnetic field. 

We explain these two phenomena within a perturbative approach known 
as mode-coupling theory. Linear stability analysis explains the instability of 
an initially flat interface to sinusoidal perturbations known as modes. In the 
initial, linear stage of pattern formation, modes grow or decay independently 
of each other. One mode, which we call the "fundamental", grows faster 
than all others. As these perturbations of the flat interface grow, they evolve 
through a weakly non-linear stage, in which modes couple with each other, 
to the strongly nonlinear late stages in which a Fourier decomposition of the 
interface shape becomes inappropriate. 

We carry out our mode-coupling expansion to third order. Linear sta- 
bility analysis explains neither interfacial symmetry breaking nor finger tip- 
splitting. At second order, we find the viscosity contrast A (defined as the 
difference between the two fluid viscosities divided by their sum) breaks the 
symmetry of the interface by enhancing growth of subharmonic perturba- 
tions to the fundamental mode. This mechanism occurs independently of 
the applied magnetic field. At third order we find a mechanism for finger 
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tip-splitting driven by mutual repulsion of elements of magnetic fluid. In the 
absence of a magnetic field, finger tips do not split in rectangular geometry 
Hele-Shaw flow. 

2 Hydrodynamics in a Hele-Shaw cell 

This section begins with a discussion of basic hydrodynamic equations gov- 
erning the motion of fluids confined within a Hele-Shaw cell, considering 
ferrofluid in particular. We present Darcy's law in the presence of a perpen- 
dicular magnetic field, and we discuss boundary conditions obeyed at the 
two-fluid interface. Since the basic equations are well established by previ- 
ous investigators, we simply review the chief assumptions and results. This 
section concludes by describing our perturbative approach. We introduce 
a Fourier decomposition of the interface shape and derive coupled, nonlin- 
ear, ordinary differential equations governing the time evolution of Fourier 
amplitudes. 

2.1 Governing equations 

Consider two semi-infinite immiscible viscous fluids, flowing in a narrow gap 
of thickness b, between two parallel plates (see figure [D). Denote the densities 
and viscosities of the lower and upper fluids, respectively as pi, T]i and p 2 , 
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r]2- Between the two fluids there exists a surface tension a. Inject fluid 1 at 
constant external flow velocity Voo = v^y at y = — oo and withdraw fluid 
2 at the same velocity at y — +00. We describe the system in a frame 
moving with velocity Voo, so that the interface may deform, but it does not 
displace from y = (dashed line in figure |l|) on the average. During the flow, 
the interface has a perturbed shape described as y = ({x,t) (solid curve in 
figure [I]) over the range < x < L in the comoving frame. 

In order to include the acceleration of gravity g, we tilt the cell so that the 
y axis lies at angle j3 from the vertical direction. To include magnetic forces, 
we apply a magnetic field Hq at right angles to the cell. By assumption, the 
upper fluid acquires magnetization M, while the lower fluid is nonmagnetic. 
We consider the limit L — > 00 to simplify calculations of magnetic forces. 

Hydrodynamics of ferrofluids differs from the usual Navier-Stokes equa- 
tions through the inclusion of a term representing magnetic force. Let M 
represent the local magnetization of the ferrofluid, and note that the force 
on M depends on the gradient of local magnetic field H. The local field 
differs from the applied field by the demagnetizing field of the polarized fer- 
rofluid. Restricting our attention to small velocity flows of viscous fluids, we 
ignore the inertial terms and write the Navier-Stokes equation for a single 
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fluid 

-rjSJ 2 u = -Vp+(M -V)H + pg (1) 

For the two dimensional geometry of a Hele-Shaw cell, the three dimen- 
sional flow u, governed by equation (|l]), may be replaced with an equivalent 
two-dimensional flow v by averaging over the z direction perpendicular to 
the plane of the Hele-Shaw cell. Imposing no-slip boundary conditions, a 
parabolic velocity profile and assuming constant magnetization parallel to 
H , one derives Darcy's law for ferrofluids in a Hele-Shaw cell 0, ||, 

b 2 f - 2M - M „ 4 ] 

The magnetic scalar potential ip is evaluated on the top plate of the cell. 
The velocity depends on a linear combination of gradients of p and ip, so 
may think of the magnetic scalar potential as part of an effective pressure. 
Equation (TJ) describes nonmagnetic fluids by simply dropping the terms 
involving magnetization. 

It is convenient to rewrite equation (0) in terms of velocity potentials 
because the velocity field v is irrotational. Since we are interested in pertur- 
bations of the velocity field around a steady flow, we write v = Voo — V <fi, 
where defines the velocity potential. Both sides of equation (fj) are rec- 
ognized as gradients of scalar fields. Integrating both sides of equation (Q) 
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yields 



1 b 2 f 2M 
rj<p= — VP V 3 + cos }> + rjVooy (3) 



after dropping an arbitrary constant of integration. 

Subtract equation (|3|) for one fluid from the same equation for the other 
fluid, then divide by the sum of the two fluids' viscosities. This yields an 
equation for the discontinuity of velocity potentials valid at the two-fluid 
interface 

.(^) + (^U^( fe -,o-^.)+,, («, 



The viscosity contrast 



A = ^i (5) 
V2 + Vi 



will play a key role in interfacial symmetry breaking. U is a characteristic 
velocity associated with driving forces, 

b 2 {p 2 - pi) g cos (3 
U= 12(„ 1+ „ 2 ) +AV - (6) 

The pressure jump across the interface, pi —p 2 depends on k, the interfa- 
cial curvature in the plane of the Hele-Shaw cell. In general this relationship 
depends upon the discontinuity of the viscous stress tensor. Under the as- 
sumption of low capillary number flow we neglect that dependence and write 
simply 

P2 ~ Pi = <tk. (7) 
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We substitute the pressure jump boundary condition (|7|) and also intro- 
duce dimensionless variables, scaling all lengths by the gap size b, and all 
velocities by the characteristic velocity a/12(rji + 7/2). The final equation of 
motion reads 



A- 



+ 



2[Uy + K - N B 1] 



\v=C- 



Ng is the dimensionless magnetic Bond number 

2M 2 b 



N 



B 



a 



and the integral 



^ — J-^oo dx' j^ x i) dy' 



1 



x — x') 2 + (y — y') 2 
1 



(9) 



(10) 



^(x-x') 2 + (y-y') 2 + l_ 
is proportional to the magnetic scalar potential. Equation (|j) governs the 
flow for a given interface shape (. 



2.2 Mode-coupling analysis 



We begin by representing the net perturbation ((x, t) in the form of a Fourier 



series 



COM) = J2Ck(t) exp(ikx), 
k 



(11) 
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where Ck(t) denotes the complex Fourier mode amplitudes. Expansion (O) 
includes a discrete (rather than continuous) set of modes k because we focus 
on the interaction of three particular modes in the subsequent discussion. 
The k = mode vanishes since we are in a comoving frame. The wavevectors 
are constrained to lie on the x axis, but can be either positive or negative. 

Now define Fourier expansions for the velocity potentials 0j, which must 
obey Laplace's equation V 2 0« = 0, the boundary conditions at y — > ±00, and 
include the discrete modes k entering the Fourier series (jTT|). The general 
velocity potentials obeying these requirements are 



'1 = 01* (*) ex P(l%) exp(iA;x), (12) 



and 



^2 = hk(t) exp(-\k\y) exp(ikx). (13) 

k^O 



In order to substitute expansions ( |T2"D and (|H|) into the equation of mo- 
tion (|8]), we need to evaluate them at the perturbed interface. For example, 
expand the lower fluid velocity potential 4>i\y=c,i evaluated at the perturbed 
interface, to third order in (. Its Fourier transform is 

Mk) = <M*) + £ mfavMtk-v + \ J2( k ') 2 Mt)( q ( k - k ^ q . (14) 

W k',q 

A similar expression for <p2\y=( can be easily obtained. These results define 
the Fourier transform of the left-hand-side of equation (pi). 



Now we must evaulate the Fourier transform of the right-hand-side of 
equation (El). The curvature in the x — y plane is M 



dx 2 



1 + 



dx 



(15) 



We expand this up to third order in £ and Fourier transform, 



k(k) = -k 2 ( k - - ( k ') 2( l[ k ~k' - q]Ck'Cq(k-k'-q- 



(16) 



The expansion to third order in powers of C, of the integral (|T0|) related 
to magnetic scalar potential, is 

1 1 



l(x) 



oo 

1 r°° 

6 J-oo 



[(x - x') 2 } 1 / 2 [{x - X') 2 + 1]V2 

1 1 



l((x') - ((x)]dx' (17) 
[C(x') - ((x)] 3 dx>. 



_[{x - x') 2 f/ 2 [{x - x') 2 + If I 2 _ 
When Fourier transformed, the integrals in ( |i~7D can be solved in terms of 
modified Bessel functions IBI 



vK > JfeT(l/2) { 1 Jo (x 2 + t 2 Y+V2 



. 00 cos kx , 
1 dx. 



(18) 



We define the functions 



J{k) = log i^\+K (\k\) + C 



(19) 



with C the Euler constant, and 



|fc|(41og2-31og3) + 2^(3^1) 

2- 

3- 



- 4#i(2|A;|) + 2# 1 (|A;|) - 
9 



(20) 



and write the Fourier transform 

l(k) = -2J(k)( k + WT(k-k'- q)Ck< q Ck-w- q . (21) 

k',q 

For nonzero k, J(k) is positive and T(k) is negative. The expansion in powers 
of ( can easily be extended to arbitrarily high order. Tsebers |J presents 
T{k) up to the fifth order term. 

To close equation (§) we need additional relations expressing the velocity 
potentials in terms of the perturbation amplitudes. To find these, consider 
the kinematic boundary condition relating the interface shape back to the 
fluid flow. The condition that the interface move according to the local fluid 
velocities is written 



dt \dx dx I r \ dy I * 



(22) 



Expand this to third order in ( and then Fourier transform. Solving for 4>ik(t) 
consistently to third order in ( yields 

(pik(t) = -jjj + J2 s 9 n ( kk ')(ki(k-k> (23) 



k' 



TT\ s 9n{k'q)C,k< q -k'C,k~ q + TT\ [ k ~ Q ~ TT J Ck'CqCk-k'-q 

and a similar expression for 4>2k(f)- The sgn function equals ±1 according to 
the sign of its argument. The overdot denotes total time derivative. 

Substitute this last expression for 4>ik(t) into equation (|TJj) for the Fourier 
transform of (f>i\ y= Q, and again keep only cubic terms in the perturbation 
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amplitude. Repeat the same procedures for fluid 2. The velocity potentials 
have now been eliminated from Darcy's law and the differential equation 
of the interface is 



= \(k)( k + A\k\ £ [1 - sgn(kk')} C k <k-k' 
+ £ \k\\q\sgn(k'q) [1 - sgn(kq)]C,k'C,q-k'C, k ~q 

k',q 

k! \k'\\k\ 
k — q — — 



(24) 



+ £*' 



2fc' 



Cfc'CgC 



qSk — k'—q 



E f^r(* -*'-(?) + ||*|(A0 a # - ^ - g] 



(k>(q(k-k> 



-q- 



Here 



X(k) = \k\[U + 2N B J(k) — k 2 ] (25) 

is the dimensionless linear growth rate multiplying the first order term in £. 
The second term in equation is second order in (, and the remaining 
terms constitute the third order contribution. 



3 Weakly nonlinear evolution 

This section analyzes the evolution of an interface under the mode-coupling 
equation ( p4[ ) derived in section [2.2| . We systematically examine terms in 



order of their strength at the onset of the instability. Thus, we begin by 
describing the first order term, which captures the well known linear insta- 
bility leading to viscous finger growth. Driving forces causing the instability 

11 



include magnetic repulsion within the ferrofluid. We move on to the second 
order term, noting the interesting coupling of a fundamental mode and its 
own subharmonic. This term is responsible for finger competition. Magnetic 
forces do not contribute to this, or any even order, term. Rather, finger com- 
petition depends upon the viscosity contrast A. We conclude our discussion 
at third order. Here, we show that finger tips may split due to coupling 
of a fundamental mode with its own harmonic. The process depends upon 
the presence of magnetic repulsion within the ferrofluid. It does not occur 
without a magnetic field. 

3.1 First order 

First order in the mode-coupling expansion reproduces conventional linear 
stability analysis. Each mode grows or decays independently of all others, 
with exponential growth rate X(k) given in equation (f2~5l). Positive values 
of X(k) make a mode unstable to growth of an initially small perturbation. 
Figures ^|- [|plot this function for three distinct cases: U — — 1, U = and 
U = +1 respectively. For each value of U we graph X(k) for three values of 
the magnetic Bond number, Nb = 0,2,4. 

In general these plots display a range of wavenumbers over which X(k) > 
0. We define two special wavenumbers: k*, the wavenumber of the fastest 
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growing mode, maximizes X(k); k c , the threshold wavenumber beyond which 
all modes are stable, is the largest wavenumber for which X(k) vanishes. 
When U — 1 and Nb = we have kf — l/VS and k c — 1. The magnetic 
field is destabilizing. As the magnetic Bond number grows, k* and k c shift to 
the right and modes of higher wavenumber become unstable. Likewise, for 
any particular mode k, the growth rate X(k) increases, causing perturbations 
to grow more rapidly. 

To analyze mechanisms of pattern selection, we will focus our attention 
on the interaction of one large amplitude perturbation, which we call the 
"fundamental" , with small amplitude perturbations of its own harmonic and 
subharmonic. We take the fundamental wavenumber kf — k* of the fastest 
growing mode. The harmonic mode kh = 2k f always lies to the right of 
the threshold wavenumber k c , so the harmonic mode is always linearly stable 
against growth. The subharmonic mode k s = kf/2 usually lies in the unstable 
regime. If present in the initial conditions it will grow, but less quickly than 
the fundamental. 

Growth of the fundamental mode creates a sinusoidal oscillation of the 
initially flat interface, forming fingers of each fluid penetrating into the region 
previously occupied by the other fluid. The interface is symmetric, with up- 
wards and downwards fingers having identical length and width. Depending 
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upon the phase of the subharmonic relative to the fundamental, either the 
upwards-pointing fingers, or the downwards pointing fingers may have their 
length modulated. The subharmonic can break the up-down symmetry of 
the growing pattern. However, within the linear stability analysis, no unique 
phase of the subharmonic is favored. Assuming the relative phase is deter- 
mined by random perturbation of the flat interface, the growing pattern will 
retain statistical up-down symmetry. For any given initial condition, sym- 
metry will be broken, but averaged over all initial conditions, symmetry will 
remain. 

Splitting of fingers is not predicted by linear stability analysis, because 
the harmonic mode is required to split fingers, and the harmonic mode is 
linearly stable. 



3.2 Second order 



Inspecting the mode coupling equation (f24|) , we note that the second order 
term does not involve magnetic field. We have previously [|T^] analyzed the 
role of the second order term in rectangular flow geometry for non-magnetic 
fluids. The results are unchanged, so we will simply recall two essential 
facts: the second order term generates finger competition dependent upon 
the viscosity contrast A; the second order term does not generate finger tip- 
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splitting. We explain these two points briefly. 

Finger competition is linked with the amplitude and phase of the sub- 
harmonic mode. Coupling of the fundamental kf to the growth of its sub- 
harmonic k s accelerates growth of the subharmonic, and selects a preferred 
phase. The selected phase varies the relative lengths of fingers of the less 
viscous fluid penetrating into the more viscous fluid. Fingers of more viscous 
fluid penetrating into the less viscous fluid tend toward equal lengths. The 
subharmonic therefore breaks the statistical up-down symmetry of the linear 
stability theory. 

Finger tip-splitting requires the harmonic mode. In the radial flow geom- 



etry [IT], second order terms drive growth of the harmonic mode kh despite 
its linear stability. These terms are absent due to the rectangular flow ge- 
ometry. As we explain in the following section, one must examine the third 
order terms to understand growth of the harmonic. 

We conclude this discussion with an explanation for the absence of mag- 
netic field effects at second order. Although the basic equation of motion (Q) 
is written in terms of forces, it is simplest to carry out the discussion in terms 
of energies. Consider the magnetic energy for a given interface shape ((x). 
The magnetic energy is unaffected by rotation of the entire experiment (Hele- 
Shaw cell and magnet) by 180° around the x axis. Because the ferrofluid is 
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paramagnetic, the magnetic energy is invariant under reversing the direction 
of the applied field. The combination of the two symmetries, rotation of the 
experiment followed by reversal of the applied field, amounts to reversing the 
sign of the interfacial displacement ((x). Since the magnetic energy cannot 
be affected by this change, it must be an even function of ((x). The magnetic 
force is given by the change in magnetic energy with respect to variation in 
interfacial shape, so it must be an odd function of ((x). 



3.3 Third order 

This section shows how the magnetic field qualitatively alters the mecha- 
nism for splitting of finger tips. We first review previous results explaining 
the general absence of tip splitting in rectangular geometry flow of ordinary 
fluids [10] . Then we describe a new mechanism for splitting finger tips in the 



presence of a magnetic field. 

We consider the influence of the fundamental and sub-harmonic modes 
on the growth of the first harmonic. Finger tip-splitting is associated with 
the magnitude and phase of the harmonic mode 2kf. It is convenient for 
the subsequent discussions to consider sine and cosine modes, rather than 
the complex modes employed in equation (0). Describing the fundamental 
as a cosine mode with positive amplitude, we only need to examine the 



16 



subharmonic and harmonic cosine modes to analyze finger competition and 
tip-splitting. Let denote the amplitude of the cosine mode of wavenumber 
k. 

Earlier papers considered finger tip-splitting for Hele-Shaw flow of non- 



magnetic fluids in the radial [i 1] and rectangular [10] geometries. Of course, 
the same results hold for ferrofluids in the absence of applied magnetic fields. 
The principal results are as follows. In the radial geometry, a term pro- 
portional to a\ drives growth of the harmonic with the phase appropriate 
to split finger tips. In the rectangular geometry, this second order term is 
missing. Instead, there is a third order driving term proportional to a>k f a\ s . 
This term is expected to be too small to split finger tips. There is also a 
reduction in the effective stability of the harmonic mode for large amplitude 
of the fundamental, but this effect cannot make the harmonic linearly un- 
stable. Consequently, finger tips do not split under normal circumstances in 
the rectangular geometry. 

Now we investigate the connection between the applied magnetic field and 
the occurrence of finger tip-splitting observed in Hele-Shaw cell experiments 
with ferrofluids. For consistency with experimental results [ffl we consider 
the case U = 1. The equation of motion for the harmonic mode (neglecting 
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terms of order 0{a\ )) is 



3 

-k h k 2 s k f [k f + 2k s }+N B k h 



(26) 



a k f a ks . 



We incorporate certain third order terms into the effective linear growth rate 



+ 



+ 



(27) 



k 2 fkh 



2 

k'ikh 



k 2 t + — A;? 



it; 



+ ^T(h) 

l -T{k s ) + -Lt(A*) 



In equations (pq) and fl27D some terms are explicitly multiplied by Nb and 
others are not. We refer to the former as "magnetic" terms, and the latter as 
"nonmagnetic" . The nonmagnetic terms reproduce the known mode coupling 
equation for nonmagnetic fluids [T0|| . 

Our mechanism for splitting finger tips focuses on X e ff- Initially, this 
quantity is close to X(kh), which is strongly negative because the harmonic 
mode is stable in the linear theory. However, A e // is increased by the presence 
of the modes kf and k s , because the coefficients multiplying their squared 
amplitudes are positive. To verify this point, recall that kf > 1 / a/3, making 
the nonmagnetic contribution manifestly positive. Also, the values of T(k) 
are negative, making the magnetic contribution manifestly positive. Since we 
consider the case in which a^ f is considerably larger than a^, the dominant 



corrections to the effective growth rate come from the terms multiplying a| . 
In the following, we concentrate our discussion on those terms. 

Nonmagnetic terms make the effective growth rate less negative but can- 
not make it go positive. The physical reason that these terms do not make 
X e ff positive can be understood by considering the contour length of the inter- 
face. Introducing the harmonic always increases the contour length, although 
the larger the amplitude of the fundamental, the smaller the increase upon 
introducing the harmonic. Multiplying the contour length by the surface 
tension yields surface energy that favors minimum contour length. Mathe- 
matically, the nonmagnetic term of order a\ f makes A e // less negative, but if 
higher orders in perturbation theory were included it would be evident that 
Ofc. cannot drive A e // positive without assistance from the magnetic terms. 

The terms in A e // that are multiplied by Nb allow A e // to eventually 
go positive, permitting growth of the harmonic. The effective growth rate 
remains negative up to a threshold value of a,k f for which A e // = 0. When a^ f 
grows beyond this threshold value, a^ h grows rapidly. The threshold value of 
a,k f should vary as the inverse square root of the magnetic Bond number, so 
tip splitting emerges sooner in strong magnetic fields. 

The harmonic mode enters spontaneously through the third order driving 
term proportional to ak f a k a i n equation (Po]). As long as \ e ff remains nega- 
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tive, this small driving force should be of little consequence. After A e // goes 
positive, however, this term can introduce a harmonic even if none is present 
in initial conditions. The existence and phase of the spontaneously generated 
harmonic depends on interplay of the fundamental and the subharmonic. 

To illustrate the occurrence of finger tip-splitting when an external mag- 
netic field is applied, we consider the interaction of modes kf and kf/2 with 
the forced modes 2k f and 3kf/2. Mode 3kf/2 behaves similarly to the sub- 
harmonic kf/2 and induces more finger competition. In figure [| we plot the 
interface evolution using the full solution to third order of equation fl2"l|). We 
examine the case in which U — 1, assuming that fluid 2 is a ferrofluid. An 
external magnetic field is applied (Ng = 1.0) normal to the cell plates. The 
initial condition is ak h = 0.145 and a>k s = —1/5 a^ r The harmonic mode is 
absent initially. Times shown are t = 0, 1, 2, 3. 

The effective harmonic growth rate X e ff starts being strongly negative 
at t — 0. It increases with time, goes through zero t « 2.16, and become 
positive. This process leads to finger tip-splitting by t = 3, as shown in 
figure [5]. Finger tip-splitting only occurs after the fundamental has grown 
sufficiently that A e // goes positive. The selected phase of the harmonic splits 
the tips of fingers whose length is variable, the less viscous fingers. 
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4 Conclusion 



Several features of the patterns formed in Saffman- Taylor experiments with 
ferroffuids can now be explained. First of all, there is the striking asymmetry 
of the interface. Since the dense upper fluid is a glycerine-based ferrofluid 
with high viscosity, and the less dense lower fluid is white spirit, we can 
understand the initial asymmetry of the interface purely on the basis of 
the viscosity contrast A, as discussed in section pT2"| . Indeed, with the field 
turned off the interface is quite asymmetric, with short and wide upwards 
fingers of the less viscous fluid and long, thin downward fingers of the more 
viscous fluid. Magnetic field effects can exaggerate an already asymmetric 
interface, but they cannot break the symmetry by themselves. It would be 
of considerable interest to perform the experiment with immiscible viscosity- 
matched fluids. 

Next, there is the splitting of finger tips which is not normally observed 
in rectangular geometry flow. Both upwards and downwards fingers are split, 
consistent with a positive value of X e ff permitting the growth of harmonics of 
any phase. The upwards fingers are more strongly split, however, consistent 
with the phase preferred by the driving force in equation (|2Bp. Both are 
predictions of our third order analysis in section pO . It would be of interest 
to examine the relationship between the onset of tip splitting and the strength 
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of applied magnetic field experimentally. 

The above features are explained by our mode coupling theory Another 
notable feature of the patterns is the nearly constant width and regular spac- 
ing of downward pointing fingers of high viscosity fluid. The constant width is 
probably the known field-dependent preferred finger width of the labyrinthine 
instability M, proportional to the plate spacing b. Given a directed set of 
thin fingers, magnetic forces will drive them towards maximal spacing, re- 
sulting in a regularly spaced array. These issues lie beyond the scope of our 
low order mode coupling approach, but may be amenable to more general 
forms of weakly nonliner analysis. 
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Figure 1: Schematic configuration of the rectangular flow geometry. The 
upper fluid is a ferrofluid. 
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Figure 2: Plot of X(k) for U = -I. 
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Figure 3: Plot of X(k) for U = 0. 
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Figure 4: Plot of X(k) for U = 1. 
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Figure 5: Plot of an evolving interface with an applied magnetic field. 
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